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Abstract 



^\ ' We explore skewed parton distributions for two-body, light-front wave functions. In order to access all kinemat- 

ical regimes, we adopt a covariant Bethe-Salpeter approach, which makes use of the underlying equation of motion 
(here the Weinberg equation) and its Green's function. Such an approach allows for the consistent treatment of the 
non-wave function vertex (but rules out the case of phenomenological wave functions derived from ad hoc potentials). 
Our investigation centers around checking internal consistency by demonstrating time-reversal invariance and conti- 
nuity between valence and non-valence regimes. We derive our expressions by assuming the effective qq potential is 
independent of the mass squared, and verify the sum rule in a non-relativistic approximation in which the potential 
is energy independent. We consider bare-coupling as well as interacting skewed parton distributions and develop 
approximations for the Green's function which preserve the general properties of these distributions. Lastly we apply 
our approach to time-like form factors and find similar expressions for the related generalized distribution amplitudes. 



1 Introduction 

Recently there has been considerable interest in the connection between hard inclusive and exclusive reactions, which 
has been, in part, due to the unifying role of skewed parton distributions (SPD's) g, Aside from being the natural 
marriage of form factors and parton distributions, SPD's appear when one tries to calculate scattering amplitudes such 
as those for Compton scattering ||, ||, the electroproduction of mesons || J?]], and di-jet production by pions on nuclei 

There is much effort underway related to the measurement of these functions |H]. Intuitively clear, but simple models 
have been used to provide first calculations (more sophisticated approaches have been pursued |l2|, [l3| Q) 
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and physical interpretations have been elucidated in the forward limit |lq ]. In this paper, we attempt to gain intuition 
about the structure of these distributions by considering two-body, light-front models as an explicitly worked example. 
Continuing with our previous work Jl6[ , we consider the SPD for a toy scalar meson consisting of two, equally massive 
scalar quarks. The features which are general in our approach are the energy denominators, which will be encountered in 
any calculation of SPD's (e.g. for a hadron of spin-i consisting of spin-i quarks). Simple two-body, wave function models 
of SPD's suffer from an inability to describe all kinematical regimes physically accessible. On the other hand, using a 
light-front Fock-space decomposition of the initial and final hadron states gives one a powerful general expression for the 
SPD's as sums of diagonal and non-diagonal overlaps ]17j encompassing all kinematical regimes. In practice, however, 
this approach is limited by our current knowledge of the higher Fock-space components of hadrons and one seems 
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constrained to use Gaussian Ansatze for these components jflj for which essential physics may be absent. Generalizing 
a rather prescient method due to Einhorn [l8| , we avoid the need to obtain higher Fock-space components explicitly by 
using the Bethe-Salpeter equation, and crossing symmetry. At this stage, we sacrifice exactness (by neglecting the mass 
squared dependence of the interaction throughout) in order to build intuition about these distributions in this approach. 
We find that one need only use the lowest Fock-space components (here the qq pair) and the equation of motion to 
access all kinematical regimes of the SPD with model wave functions. We note that recently a similar analysis, although 
tentative, appeared for spin-i quarks in the context of SPD's in [ ^9[ . 

We begin in section |^ by writing the covariant Bethe-Salpeter equation for our meson and extract from this equa- 
tion the light front wave function and its equation of motion (the Weinberg equation). By way of orthogonality and 
completeness, the Green's function for the Weinberg equation is then defined. In section || we apply this analysis to the 
covariant, leading-twist diagram for deeply virtual Compton scattering (DVCS). We carefully demonstrate that in the 
deeply virtual limit, SPD's can be calculated from the light-front time-ordered triangle diagrams for the form factor, 
making clear the close connection between the sum rule and light-front time-ordered perturbation theory. These dia- 
grams are evaluated in section ^ in terms of qq light front wave functions where the Z-graph is dealt with by crossing the 
interaction. The resulting bare-coupling SPD's are derived and their simplicity is used to discuss continuity, x-symmetry 
and time reversal invariance properties. The demonstration of time reversal invariance for the bare-coupling SPD's is 
started in section |] and completed in Appendix A. 

The bare-coupling SPD's are only (partially) valid at high \t\, which is unfortunately not likely accessible in experi- 
ments. Moreover, the effects of higher Fock components in this approach are contained in the photon vertex function. 
Thus inclusion of interactions at the photon vertex is required and carried out in section ||, where we derive the photon 
vertex function in terms of the interacting theory Green's function. Finally the general form of the SPD's is deduced in 
section ^. These distributions satisfy the necessary symmetry properties and are continuous due to the behavior of the 
Green's function near the end points. As calculation of the Green's function is likely burdensome, we discuss in section 
possible approximations to the Green's function which maintain the symmetry and continuity of the distributions. 
We explore first the Born approximation valid at high \t\ (which does not result in the bare-coupling SPD's). Next we 
examine the non-relativistic scheme and finally the closure approximation. 

The stage is set to apply the non-relativistic approximation scheme to the Wick-Cutkosky model in the weak-binding 
limit (section ||) . We use this example to illustrate how to cross the interaction for partons moving backwards through 
time and thereby explicate that in general an ad hoc potential (which might generate some phenomenological hadronic 
wave function) lacks the field theoretic nature to accommodate such crossing. We then verify the sum rule for the 
form factor in the non-relativistic scheme, which is accomplished by using only the valence contribution to the SPD. 
Appendix B makes explicit the recovery of the Drcll-Yan-West formula for the form factor in the limit of zero skewness. 
Furthermore, we comment there on the general consequences and limitations of approximating the interaction to be 
independent of the mass squared. 

The extension to time-like form factors is carried out in Appendix C. We proceed analogously with our consideration 
of SPD's by deriving the sum rule for time like form factors in terms of generalized distribution amplitudes. Then we 
calculate the generalized distribution amplitudes in terms of the relevant light-front time-ordered graphs. Lastly we 
conclude briefly in section pi 



2 Light-front reduction of the Bethe-Salpeter equation 

We start by writing down the equation for the meson vertex function T. It satisfies a simple Bethe-Salpeter equation, 
first schematically as operators (see Figure [j]): 

r = -VG T, (1) 

where G is the full, renormalized, two-particle Green's function (which can be expressed as a product of two, renor- 
malized, single-particle propagators). Here V is the interaction potential, which is the negative of the irreducible 
two-to-two scattering kernel. Let us denote the meson four-momentum by R, the quark's by p, and the physical masses 
by R 2 = M 2 ,p 2 = m 2 . Since our quarks are scalar particles, the renormalized, single particle propagator has a Klein- 



2 



V 



Figure 1: Diagrammatic representation of the Bethe-Salpeter equation. The blob represents the vertex function T. 



Gordon form: 

= (fc 2 - to 2 )[1 + (fc 2 - m 2 )/(fc 2 - to 2 )] + it ^ 

where the residue is one near the physical mass pole and the function /(fc 2 — m 2 ) characterizes the renormalized, one- 
particle irreducible self interactions. It is a matter of later notational convenience that we write / = /(fc 2 — to 2 ). We 
shall assume there are no poles (besides at the physical mass) in the propagator. This is consistent with model studies 
p2f . In the momentum representation, the Bethe-Salpeter equation then appears as 

r( P) R) = -iJ t^* 3 ^ WM> ( 3 ) 

in which we have defined the Bethe-Salpeter wave function ^bs as 

^ BS (k,R) = G (k)T(k,R)G (R-k), (4) 

The Bethe-Salpeter wave function is covariant, thus in order to perform the light front reduction we must project 
out the x + — hypersurface where the initial conditions, wave functions, etc. of our system are defined [^0). (We define 
the plus and minus components of any vector by x = (x° ± x 3 )/y/2.) This projection onto the initial surface gives 
us the light front wave function, vis ^ lf = / dk~^BS- The remarkable property of light front wave functions is that 
they depend not on the total momentum of the system but on the relative momenta of the constituents |2l| . Since 
the constituents of our meson are two equally massive quarks, the relative momenta take a simple form. Writing this 
dependence out explicitly, we see 

* LF (k, R) = R+ Jl_ x M *> k ^ el = k± -^ R± ) . (5) 

with k + / R + = x and the pre-factors chosen with malice aforethought. In order to perform the x + = projection on 
the Bethe-Salpeter wave function, we must know something about the analytic structure of the vertex function. We 
would like to perform the k~ integration simply by assuming there are neither fc~ poles, nor branch cuts in the vertex 
function. Looking at the effective qq potential in Eq. (||), shows that the p~ dependence of the interaction V(k,p) gives 
rise to the p~ poles of the vertex function T(p, R). Thus taking the effective qq potential V(k,p) to be independent of 
the minus momenta eliminates the vertex function's dependence on the light-front energy and allows us to easily perform 
the light-front projection. The effective two-body interaction in 1 + 1 QCD, for example, is independent of the minus 
momentum. Though certainly this independence is far from general. Calculating V(k,p) from light-front, time-ordered 
graphs, however, requires the on shell condition which removes any functional dependence on the minus momenta from 
V and hence T. Now carrying out the k~ integration on the Bethe-Salpeter wave function, we find 

V>(z,k\. ol ) = R+x 0-_- x ) f dk-^ BS {k,R) 
m J 

= A w (x,k ± rel |M 2 )r( 2 ;,k i 1 , 1 ), (6) 
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where we have defined 



„ , ,± W^relPg) 

A w (2;,k re i|M ) = -, (7) 

1 +xD w - 1 (^k^ re i iM^/^Dw-^k^ej \M 2 )j 



with the help of the abbreviation D\y for the Weinberg propagator 

I 

M 2 



D w (x, k x rel \M 2 ) = - ^_ k± ^ . (8) 

x(l— x) 

The projection carried out in equation ^ can also be utilized to establish the symmetry relation 

x /(xDw^.k 1 |Af 2 )) = (1 - x)/((l - i)Dw(l - s.k 1 \M 2 )), (9) 

required for a system with equally massive quarks. 

Notice we have implicitly altered the functional dependence from T = T(k, R) to T — T(x, k rei) - As a result of the 
integration in Eq. (^|), x is restricted: x G (0, 1). Thus the light-front wave function vanishes for x outside this range. 
Eq. (^) gives the relation between the vertex function and the light front wave function. This relation will be exploited 
to evaluate Feynman graphs in terms of wave functions below — provided, of course, x € (0, 1). 

It remains to find the equation of motion satisfied by the wave function ip( x i k r e i). We accomplish this by performing 
the k~ integration in the Bethe-Salpeter equation (||). 

\ yr - 1 (x,k L tfil \M 2 )^{x,)c L tA ) = -i / £zi*Bs(p,R)V(p,k) 



(2tt) 

V(y,p re i;a;,k ro i)V'(y,P re i)> ( 10 ) 



2(2tt) 3 i/(1 - V) 



with y = p + /i? + , and p r el = p —yR . The above equation can be written as the Weinberg equation 23 by 



relocating the renormalized self interactions into the definition of the potential. To do so, we write the above equation 

* k i2 +m 2 

more suggestively as an eigenvalue equation: (M 2 — H)ip n — 0, where H = H Q + V with H Q as the kinetic term 



x(l — x) 

and V as the potential, whose action on ip is spelled out by 



V^(x,k A 

with 



dydp A 



2(27r) 3 t/(l - v) 



V(x, k x ; y, p- 1 ) - 2(2nfS(x - y)S 2 (k ± - p x ) A(x, k x ) 



iKi/.P" 1 ), (11) 



A(x,k- L ) =xD w - 2 (a;,k- L \M 2 )fUTh N ~\x,* L W 2 ))- (12) 



At this level, the action of the potential is still exact, however, the self interactions as well as the potential function 
depend on M 2 . If we neglect this dependence (as is possible in 1 + 1 QCD, or to leading order in a weakly bound system 
where M 2 « 4m 2 ) then the Weinberg equation becomes a standard eigenvalue equation. Under this assumption, the 
eigenfunctions {ipn} satisfy orthonormality and completeness, which appear in our normalization convention as 

dxd k 1 - 

^ n (x,k ± )^(x,k ± ) = S nm , (13) 



2(2tt) 3 x(1 - x) 

Y, Ms, ^Xiy, P X ) = 2(2tt) 3 .x(1 - x)5(x - y)6 2 (k ± - p x ). (14) 
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Figure 2: The leading DVCS diagram 



The Green's function for the Weinberg equation 

Gfek^.p^IZ ^^f' (15) 

n n 

thus satisfies the equation 

(R 2 - H^jG(x, k ± ;y, p ± \R 2 ) = 2(2tt) 3 2;(1 - x)S(x - y)6 2 (h ± - p^). (16) 

Clearly approximating away the M 2 dependence of the interaction has limitations. We spell out some of the relevant 
consequences in Appendix B. 

3 Deeply virtual Compton scattering 

Skewed parton distributions show up, for example, as generalized form factors in deeply virtual Compton scattering. 
We will introduce the SPD's in this way. In our toy model, we imagine a virtual photon interacting with our meson 
(labeled as tt) but leaving it intact P 2 = P' 2 = M 2 : 

1 *(q)+7r(P)^ 1 (q')+7r(P'). 

This process is certainly not imaginable in any realistic experimental setup, however, it is related by crossing symmetry 
to the physical reaction 7* + 7 — > tt + tt which was investigated in f24j] . The relation between SPD's and the generalized 
distribution amplitudes (GDA's) — which characterize two pion production — was pursued in [P5[ . In Appendix C, we 
illustrate this connection in our light-front approach by calculating our model's GDA via the time-like form factor. 

Returning to Compton scattering, the momentum transfer suffered by the meson we denote by A M = P'' 1 — P p (and 
its invariant t = A 2 ), where the \i labels (—,+,1,2). The leading-twist (handbag) diagram for this process is depicted 
in Figure |^ where we give the parameterization of the intermediate momenta. We follow the standard simplification []l3[ 
and choose a frame where q + = [Efl, so that q 2 = — q ±2 = — Q 2 . In the deeply virtual limit, Q 2 is taken to be large 
while the ratio -£p-^ = £ is finite. This defines C as the Bjorken variable for deeply virtual Compton scattering. In this 
section, we work in the asymmetrical frame where the initial meson's transverse momentum vanishes and consequently 
q~ = Q 2 /2(P + . Lastly, we define X as the ratio of the active quark's plus-momentum to that of the initial meson: 
A = k+/P+. 

At leading twist, we can ignore interactions at the quark-photon vertex as well as quark self-interactions between 
photon absorbtion and emission. Evaluating the diagram shown in the figure thus yields 

M%" =i J H (2k \t+ ) q ( p k -l^+te G °( fc ) F ( fc ' P ) G °( fc + A ) F *( fc + A ' P + A ^ P ~ fc )' < 17 ) 
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p p> 

Figure 3: The triangle diagram for the electromagnetic form factor 



where A signifies that this contribution is that of Figure We also need to calculate the crossed diagram Mg" to get 
the covariant amplitude . 

The denominator in ( p~T| ) which depends on the momentum q can be dramatically simplified in the deeply virtual 
limit. Furthermore there is a cancellation in this limit from the numerator of M^ + : 

+q- 1 1 

(k + q) 2 -m 2 +ie 2P+X-( + ie V ' 

There is an analogous cancellation in the deeply virtual limit when we evaluate the crossed diagram + . 

(k + A - q) 2 - to 2 + ie 2V+ X -ie 



(19) 



The remaining terms in the integral are the same for each diagram. Thus in the deeply virtual limit, we can decompose 
the Compton amplitude as 

M-+w- / 1 (2X-C).F(X,C,t)(-l-+ \.. )dX. (20) 
J \X — ie X-Q + ieJ 

Removing the dk + integral from the expressions for and MZ + reveals 

f dk^dk 1 - 

^(X,C,t) = -t J (2 ^ 4 G (fc)r(fc,P)G (fc + A)r*(fc + A,P + A)G (P-fc). (21) 

This is just the integrand of the covariant triangle diagram which gives the electromagnetic form factor. In this way we 
rediscover the sum rule of Ji 

—AT(XX,t)dX = F(t), (22) 

where the factor (2 — C) _1 appears because we normalize to the average meson momentum P M = (P' M + P M )/2. The 
kinematics gives the relation P + = P + (l — C/2). The other factor, 2X — (, removes the derivative coupling at the photon 
vertex from the SPD. The (^-independence of the form factor is due to Lorentz covariance. Once we have expressed 
JF(X,C)t) in terms of light-front wave functions, we will be concerned with verifying Eq. (E3). 



4 Triangle diagram and the Z-graph 

We have reduced determining the SPD to a calculation of the covariant triangle diagram for our meson's electromagnetic 
form factor. In the previous section, we employed a special reference frame in which the virtual photon is space-like 
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(q + = 0). Thus in calculating the triangle graph from Eq. (pi]), we no longer have the freedom remaining to choose 
A + = which is often done in calculating form factors. We must proceed on a course parallel to Sawicki J27J, who 
considered such general expressions for form factors earlier. 

The k~ integral appearing in Eq. ( pl| ) is convergent which enables evaluation by choosing to close the contour in 
either the upper- or lower-half complex plane. Without any loss of generality, we choose A + < 0. The poles of our 
integrand are located at 

K a - -A H 2(k++A+) 2(fe++A+) 

K =4^-1^ (23) 

V — r 2(P+ -k+) h 2(P+ -k+)- 

When k + < 0, all poles lie in the upper-half plane. On the other hand, for k + > P + all poles lie in the lower-half 
plane. Then for both of these cases, we merely choose the contour which avoids enclosing any poles and, by Cauchy, the 
integral vanishes. To evaluate the k~ integral for < k + < — A + and — A + < k + < P + , we must enclose at least one 
pole when closing the contour. 

First let us look at the case — A + < k + < P + . This is the familiar region which survives when A + = (of course 
now A + = — C-P + )- Closing the contour in the upper-half plane we enclose only the pole at k~ and the integral is 
+2iri Res(/c~). Dropping pre-factors for the moment, we find (see Figure^): 

T(X, k- 1 -X P x ) r*(y, kr 1 -Y P x +(1 - Y)A ± )) 

(K - k-)[l + 2k+(fc c - - fc-)/(2k+(fc c ~ - k-))} (K - ka)[l + 2(k+ +A+)(fc c - - ka)f(2(k+ +A+)(fc c 7 - *-))] 

= 4(1 - C) P + 2 i>(X, k x -X P ± )ip*(Y 1 k- 1 -Y P 1 - +(1 - Y)A ± ), 

where Y = > 0. Now we merely shift the k x integral and work in the symmetrical frame [p8| in which P^ = 0. 



i-C 

Thus for X > £ we have 



^(x,C,t)= VM * ^ M / 4^^(xMW(yM +(1 - x)^), (24) 



X(1-X)(X -C) J 2(2tt) 

where the average of the plus momentum arguments between the wave functions is x = (X + Y)/2. We have been 
employing ratios of plus momenta with respect to the initial meson. To instructively counter this malady, we change 
variables from Radyushkin's || to those of Ji To this end we define a new ratio of plus momenta with respect to 
the average meson momentum x = k + / P + , inverted this reads X — Furthermore C = The new symmetrical 
SPD (for x > £) is just F(x, £,t) = (1 + ^)!F{X{x, £), C(0i an d absorbs the Jacobian factor lurking in the sum rule 
(p2|), vis. r_ iF(x,^,t) efe = F(i). Carrying out this conversion on ( p4] ) 



(25) 



we find our earlier result (r| |2^] for x > £. 

Now when < k + < — A + , we close the contour in the lower-half plane and the integral is — 2m Rcs(fc^). Once again 
temporarily ignoring the pre-factors, the residue gives (see Figure |J) a wave function from the initial meson's vertex: 

Res(fc~) oc ^(XM +XA 1 /2) Aw^k 1 +(1 - Y/2)A X \M 2 )T*(Y, k 1 - +(1 - Y/2)A X ). 

The final meson's vertex, however, represents meson production off the quark line in this time ordering and thus is a 
non-wave function vertex. Said another way, since Y = < 0, Eq. (^]) does not hold.[|] The vertex function, on the 
other hand, is defined for Y < since the physics of a non-wave function vertex is related by crossing symmetry to that 

1 We could define a non- valence wave function from the non- wave function vertex but Eq. ^6| shows how this non- valence wave function 
would be related by crossing to the valence one. 
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Figure 4: Time-ordered triangle: the Z-graph 




Figure 5: How the interaction crossing symmetry fixed up the non-wave function vertex. 



of the wave function vertex. Using the Bethe-Salpeter equation (^|) for the vertex function, we can perform the crossing 
by inserting the interaction 



(26) 



which is well defined since only the interaction contains the negative momentum fraction. The interaction has enabled 
us to cross lines and recover a wave function at the final meson's vertex, see Figure |[ We will comment on dealing with 
the crossed interaction in light front field theory (using the Wick-Cutkosky model as an example) below in section |[ 
Having crossed lines, we arrive at 



^(x,£t) = 



-c 



dk^dn^dZ 



X(1-X)(X-C) J [2(2tt) 3 ] 2 Z(1 - Z) 

x A w (a, k 1 - +a A 



i>(X,k x +IA 1 /2) 

t)v(z, K x ; Y, k x +(1 - F/2)A ± ) r(Z, k x ), (27) 



with a = X/(. 

Without interactions at the photon vertex, Aw is an approximation to the qq component of the photon's wave 
function. When £ — > 0, this approximation does not cause any harm since we recover the Drell-Yan-West formula 
(see Appendix B). Additionally, the above expressions guarantee that the SPD is continuous at the crossover (X = Q. 
This is because for a physically reasonable potential (one which vanishes at least linearly at the end points), Eq. (10) 
mandates the wave function vanishes quadratically (as pointed out in (3^]). Thus Eq. (HJ) vanishes as X — * £ due 
to the wave function and Eq. (|27j ) vanishes in the same limit due to the interaction (in both cases the offending 
X — C in the denominator is canceled and what remains vanishes linearly). Continuity is required by the factorization 
theorem for the hard DVCS subprocess. Additionally, due to the reality of the leading-order Bethe-Heilter amplitude, 
Im(M'"') oc F(C, C, t) is physically observable [^0j which is another reason continuity is essential. 

Treating the photon vertex as bare, however, requires high |f|, which is not only an extraneous assumption but also 
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not experimentally realizable.^ Furthermore for the sake of the sum rule, we must make sure the effects of all the Fock 
space components are present in this non-perturbative scheme. To do so, we must include interactions at the photon 
vertex. 

While the bare-coupling SPD's derived above aren't of general applicability, they are simple enough to demonstrate 
the required symmetry properties of SPD's, and hence shall not be immediately discarded. To form the SPD for our 
flavor singlet meson in the symmetrical variables, we use 0] 



(1 + C)- 1 F(x, & t) = ^(X, C, t) 0[X(1 - X)} + T(( - X, C, t))6[(C -X)(X + 1-Q], 



(28) 



with which we can test time reversal invariance which requires invariance under £ — > — £. Firstly, however, it is 
instructive to investigate F(x, £, t) for x < — £ which is !F(C — X, £, t) for £ — 1 < X < 0. Straightforward algebra shows 
F(x,t;,t) = F(\x\,£,t) for x < — £. Furthermore, the form of equation (25) shows F(x,—£,t) = F(x, £,t) in accordance 
with time reversal invariance. These observations show we are on the right track. 

In the central region — £ < x < £ the £ symmetry is obscured. Writing out the consequence of equation (28) in the 
central region, the x — > —x symmetry is obvious 



F(x,£,t) 



-2£(l + 



dk^dK^dZ 



(x 2 - £ 2 )(1 - x) J 4(2tt)6Z(1 - Z) 



■1>(x, k x +1AV2) A w (I + k - 



2C 



1 x . 
•2 + 2^ 



x V 



(z, k^; Y, k^ +(1 - F/2)A ± ) $*(Z, K X ) + I 



x -> -a; ^ (29) 



and follows directly from Eq. fl28|). Note that we now use X and Y as replacements for yi| and j5f> respectively. 
When we separately calculate the expression for the time-reversed process in Figure ^, we arrive back at the complex 
conjugate of (|2^). On the other hand, substituting £ — > — £ in Eq. (^9|) in order to verify time-reversal invariance is 
fraught with difficulty. Evoking this transformation sends X — ► Y and Y — * X, and formally the above expression is 
zero since ip(Y < 0, . . . ) = 0. To investigate time-reversal invariance via £ symmetry, we must do so before appealing to 
crossing symmetry, since in the reverse course of events the wave function and non-wave function vertices are switched. 
Furthermore one must reanalyze the locations of the poles, etc. Indeed after much algebraic mancevering, one recovers 
Eq. (^9|) for A — > — A (the details are relegated to Appendix A). At this point, we have demonstrated the x symmetry as 
well as the time- reversal invariance of the bare-coupling SPD's. The x symmetry will obviously extend once we include 
interactions due to the form of Eq. ^8|. On the other hand, showing the £ symmetry will require a slight reworking, 
though the scheme set up in the appendix is easily extended. 



5 Interactions at the photon vertex 

For general applicability at all values of t, we saw above that interactions at the photon vertex must be included and so 



we derive the T- matrix for qq scattering by adapting the method in 31 



The T-matrix satisfies a four-dimensional Lippmann-Schwinger equation. We are interested in its three dimensional, 
light-front reduction which gives us the equation 

T(x, k ± ;y, p 1 - \R 2 ) = (y, p ± ;R 2 \ f \x, k x ; R 2 ) 

/dzda 1 ^ 
2(2Trfz{l-z) V{x ' k±; Z ' q±) Aw(z ' q± |i?2) T(Z ' q±; 2/1 P± |jR2) ' (30) 

where R 2 is the invariant mass of the system and we have lazily dropped the relative label from all transverse momenta. 
Now the function <f> defined by 

cf>{x, k x ; y, p x \R 2 ) = A w (x, k x \R 2 ) T(x, k x ; y, p x \R 2 ) (31) 



2 The bare-coupling vertex is actually of questionable validity even at high \t\. This is because the Born approximation to SPD's does not 
result in the bare-coupling SPD's. See section ^. 
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satisfies the inhomogeneous equation 

(ii 2 - H)<f>(x,-k s -,v,p ± \&) = -V(x,k ± ;y,p x ). (32) 
Using the properties of the Green's function ([l5|), we can construct the solution of this equation, namely 

/dzda 1 - 

and in turn write a familiar new expression for T 
T(x, h ± ;y, |i? 2 ) - -V(x, k x ; y, p x ) 

Lastly we can render the above expression for T more useful by utilizing the equation of motion for the Green's 
function Eq. ( |l6| ) to rewrite equation ^ as 

<f>(x, k x ; p x |i? 2 ) = 2(2n) 3 S(x - y)8 2 (k ± - p x ) - G(x, k x ; y, p x \R 2 ) A w -\y, p x |i? 2 ). (35) 

Using this information, we arrive at the final form for the T-matrix which will be relevant for quark-quark scattering at 
the photon vertex 

Tfok^y.p 1 - \R 2 ) = Aff-'^k 1 |i? 2 ) 

x ^2(2^) 3 x(l - s)«(a: - ^(k 1 - p x ) - k x ; y, p x |i? 2 ) Aw" 1 ^, p X |i? 2 )^) . (36) 

To add the necessary quark-quark interactions, we replace the bare photon coupling with the photon vertex 
function T^, which can be written in terms of the T-matrix: = + D^G Q T. To calculate the vertex function, 
we look at the plus component and perform the light-front reduction. Let the initial quark have momentum p and the 
scattered quark p + A. Define W — p + / P + , where P + is some external plus momentum, and keep t = A 2 as above. 
We thus have 

r+(p,-A)/iP+=(2W-C)-C / ^Qnyxtt-X)^ (2X ~ Aw(g ' k± m{a ' k±]U1 ' P±rcl ^ (37) 

where p ± Tcl = p x +uA ± , X = k+ /P+ and A+ = -CP + , cr = X/( and w = W/C- Notice that as a result of the 
light-front reduction, X is restricted to (0, £). Now we merely insert the T-matrix derived above Eq. ^to discover 

r+(p,-A)/P+ = -iCA w - 1 (a; ) p ± rel |t) J ^^x^-X)^ {2X C)G(g ' k± ; P±rel (38) 

In particular, we should note that the additive bare-coupling term canceled. The above equation tells us how to add qq 
interactions within the light-front reduced Bethe-Salpeter approach. 

6 SPD's 

In the non-perturbative Bethe-Salpeter approach to SPD's, we require quark-quark interactions at the photon vertex. 
Above we have derived an expression to include them in this approach. Thus we now proceed to derive the SPD's using 
everything developed so far. 

Following the above, we merely insert the photon vertex function into equation [2l], now relabeling k as p since there 
are two loop integrals (the k loop integral is contained in T + and X = k + / P + is fixed, while W = p + / P + is integrated 
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over) . There is no new p~ dependence introduced and consequently the integration proceeds exactly as in section ^ and 
hence there are two contributions to assess depending on W. 

For < W < C, we can get the contribution from Eq. ^7] by inserting the form of T + Eq. |38| — being careful, of 
course, to remove the constants, derivative coupling, as well as the integral over X. This gives a contribution to the 
SPD for < X < C 



C < t] = -Xjhn J [2(2, W wl-W)(w UU- Z) m P +WA /2) 

x G(u, p x +uA x ;a, k x \t)V(Z, k x ; Y, p x +(1 - y/2)A x )^*(Z, /s x ). (39) 

where now we have used y = ^3 . 

On the other hand, appealing to Eq. ^4] for £ < < 1 and inserting r + leads us to G(u> > 1, . . . ). In order to 
interpret the negative momentum fraction (here 1 — u>), we must first use the equation of motion for the Green's function 
© & 

A w -\uj, p x |t)G(w, p X ; cr, k x |t) = 2(2tt) 3 A(1 - a)S{X - W)S 2 (k 1 - - p x ) 

/c\7,f] K 1 - 
2(2^(1 -Z) V ^ ^ k± I*)' ^ 

after which the problematic momentum fraction has been relocated solely in the interaction. As before, the definition 
of V can be extended in our light-front field theory by crossing. In the first term above, the delta function mandates 
£ < X < 1. While in the second term, we have < X < £ as originally required by equation 38. Indeed the effect of 
the first term is to regenerate the bare photon coupling, which gives a contribution ^(X, £, t) identical to Eq. 24 above 
0. The second term in Eq. |^ gives rise to the final contribution 

C(l-C) [ dk ± dWdp ± dZdK ± 8l(l-W)(W-0] /mr , WA± 



^ X ^ = xf-X-)J mnfiwd -W)(W- Ozfl - Z) ^ W > P +T ^ A /2) 

x V(Z, k x ;lj, p x +^A x )G(Z, k x ; <t, k x |t)V* (y, p x +(1 - F/2)A X ) . (41) 

Combining these results, we can finally write the skewed parton distribution as 

^(X, C, t) = T a (X, C, t)8[[X - C)(l - X)] + (> b (X, C, t) + W C, tj)6[X(C - X)}, (42) 

which is the 3 + 1-dimensional generalization of the SPD's found in JT^]. We should note: given the form of the Green's 
function (p"5|), the continuity of the SPD above is maintained at X — £ due, this time, solely to the nature of the wave 
function at the end points. The graphical representation of the SPD for X < ( is shown in Figure ^ 



7 Approximating the Green's function 

With Eq. ^2], we have deduced the form of the SPD satisfying the necessary symmetry properties as well as the continuity 
condition at the crossover. It remains to check whether the sum rule, Eq. |2^ is satisfied, however, evaluating an eight- 
dimensional integral for the SPD (or nine-dimensional for the sum rule) is not exactly the best way to proceed. In 
this section, we investigate possible approximations for the Green's function which are still consistent with the general 
properties of the SPD. Of course, there is the obvious simplification if t is near a bound state pole. But in general, one 
seeks to avoid solving for the entire spectrum as well as the eigenfunctions of a given light-front potential. Additionally, 
there is the Born approximation for high \t\ which for JF 6 (Eq. ^) results in the bare-coupling SPD encountered above. 
There is also the contribution from T c (Eq. |l]) at high \t\ which is absent from the bare-coupling SPD. 

3 The simplicity of the SPD for X > f in this approach is related to the expression for the form factor. We make this connection explicit 
in Appendix B 
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Figure 6: Graphical representation of contributions to the SPD for X < £: T\> + T c . The shaded box represents the 
crossed interaction V(l — > 3), and we implicitly include a sum over all meson states for the internal double line (thus 
effectively the double line represents the four-point Green's function). 

7.1 Non-relativistic approximation 

In the non-relativistic approximation, the integral equation for the T-matrix becomes algebraic. Thus one can obtain 
an approximation of the Green's function by way of the T-matrix. Let us return to Eq. (^) for the T-matrix. In a 
non-relativistic scenario (m — > oo), we may write A\y — ► Dw assuming the self interactions become negligible in this 
limit, i.e. m 2 /(m 2 ) — > 0. And thus 

/dzd q 
2(2tt)3 Z (1 - z) Dw(z ' q± |i?2)T(z ' q±; 2/1 P± |i?2) 
= -V(x, k x ; 1/2, 0)1(2/, p x \R 2 ). (43) 

The resulting unknown integral I(y, p |i? 2 ) can be solved for algebraically. The result leaves us with the above approx- 
imation to the T-matrix with 

I(y, p X \R 2 ) = - -i „ p2 . / , D w (z, |i? 2 )F(z, q x ; y, p x ), (44) 



1-C(m 2 ,i? 2 ) J 2(2tt) 3 z(1 - z) 

where 



r """'- ^ = I 2(2M1- z) BW(Z ' q± |fl2)F(Z ' 1/2 ' 0) ' (45) 



Further simplification of Eq. 44 in the non-relativistic limit is not possible, since doing so results in a logarithmically 
divergent expression. 

Inserting the approximate T-matrix into equation |36|, we can in turn produce the non-relativistic approximation to 
the Green's function 

G^^k 1 ;!,,? 1 \R 2 ) « 47r 3 <%- l/2)<5 2 (p ± )/(x, k x |i? 2 ) D w (x, k x |i? 2 ), (46) 

where is it understood that this expression is to be used where the variables y and p^ are integrated over. Notice that 
as x — ► 1, the propagator vanishes which is only enough to yield a finite SPD at the crossover. The additional term 
k \R 2 ) depends upon the interaction. Using general properties of V from the Weinberg equation and electromag- 
netic form factor, one can show I(x = 1, . . .) = |33] |. (This is just a way to motivate the physicality of requiring V to 
vanish linearly at the end points.) This feature applied to SPD's maintains continuity at the crossover between valence 
and non- valence regimes (X — £). 

In the non-relativistic approximation scheme (NR), we arrive at the distributions 



X NR + ^ NR = 



= (C h ((,t) + C c ((,t)) I ( J^? c _ x) Jfa k x \t) Dw(a, k 1 \t) (47) 



2(2tt) 3 X(C - X) 
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where we have defined 



C(l-C) [ dWd^6[{W-0{l-W)]^ r x 



Cc(c ' i) = %#7 #(i - - c) ^ p / 2 ) ni/2, o; », p ± (y, p-+(i- y/2)A ) , 

(48) 

and we have suggestively written 

« r = °) = /3p$^)^ k± )- < 49 » 

7.2 Closure approximation 

The non-relativistic approximation to the Green's function is an important limiting case but not of general applicability. 
Here we use the closure approximation which can be applied to the SPD of any system by calculating two parameters. 
Into the energy denominator of the Green's function, let us introduce the parameter (/i 2 ) 

1 1 



R? - Ml R? - ( M 2) + (ft) -Ml 

1 </?) -Ml 



R 2 - (v 2 ) {R 2 - (n 2 )) 2 



(50) 



Now requiring the contribution to the SPD from the second term in the Green's function expansion to vanish puts a 
restriction on (fi 2 ). Of course there are really two parameters (fi 2 } b and (n 2 ) c which could be calculated from having the 
resulting integrals in Eqs. (39, 4l|) vanish. Exact calculation of the parameters, however, requires the exact Hamiltonian. 
But in principle, the expansion is well defined. 

Keeping only the first term in equation [5^ and using the completeness relation (|lj) we arrive at the distributions in 
the closure approximation (CL) 

jrCL = C f d^dZd* 1 - -^( X ,k ± )V(Z, K ± ;Y 1 k ± -(l-x)A ± U*(Z, K ± ) 

t-(fj?) b J [2(2tt) 3 ] 2 x(i - x)(c- x)z(i- zy v ' ; v ' ' ' v ' y ' ' 

(51) 

where above Y = -jE^> % = (X + Y)/2 and y — . The finite- ness of the above expressions depends solely on 
the limiting behavior of the crossed potentials. Given a linearly vanishing potential at the end points, the closure 
approximated SPD is finite at the crossover, but not zero and hence the SPD is discontinuous. This is not surprising 
since the closure approximation to the Green's function is highly similar to a Weinberg propagator but without the 
necessary end-point behavior (which forced the bare-coupling SPD's to vanish). It would be possible in a given theory 
to use the end-point behavior of the crossed interaction along with the continuity condition to derive a constraint between 
(/x 2 ) b and (/i 2 ) c . Consequently the sum rule could then be exploited to evaluate the remaining parameter. This will not 
be pursued here, but will likely be the source of continued investigation. 



8 Wick-Cutkosky model wave functions 

One cannot hope to use wave functions which are not derived from field theory to verify the sum rule in this approach 
since we would be at a loss to properly extend the potential's definition. Thus model wave functions derived from ad 
hoc potentials (such as those found in our earlier work or the constituent quark model wave functions) cease to be 
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Figure 7: Diagrammatic representation of the one boson exchange potential: V = V a + Vb. 

useful. Here, we use the Wick-Cutkosky model and begin with a weak binding solution in the non-relativistic limit |$3]. 
Certainly this choice is motivated didactically since it results in clean analytic expressions with wave functions about 
which we have great intuition. Indeed the point of this paper is to demonstrate proof of concept (via time-reversal 
invariance and the sum rule) rather than trying to exploit features of SPD's to investigate model properties as was 
attempted in [pi) . 

Applying the Wick-Cutkosky model to our problem is simple, we merely use an interaction proportional to y/aq 2 (f) in 
our scalar quark's Lagrangian. Here, is a massless scalar field mediating the interaction. For a suitably weak coupling 
a, the potential can be well approximated by single <f> exchange. Summing the two light-front, time-ordered, one boson 
exchange diagrams (seen in figure |^) gives rise to our potential. Let the total momentum of the system be P and define 
z — p + / P + and y = k + / P + . Using the time-ordered rules (see e.g. we arrive at 



v = 2 fl(z-y) i ( , 9] 

a 9 2(p+ - k+) P--p--(p- k)- - (P - k)- { ' 

v = 2 °(y- z ) i fttl 

6 9 2(k+-p+)P--fc--(fc-p)--(P-p)-' [ ! 

where V a corresponds to the first diagram in figure [?] and H the second. Here g 2 = 16-Km 2 a. The sum of these two 
terms can be combined into 



with 



V{z,\^;y, V X ) = -9VK>\ (54) 
^ = k±2 S^y) _ 2 k ± .l S^z) + ^ 

z(i-z) 2/(1-2/) 

where M 2 is the total invariant mass squared, and we have again lazily dropped the relative label from all transverse 
momenta. Additionally we have defined the replacements 



2(1 — z) y(l — y) m 

z i l -y), z >v 

y{^-z), y>z 



S[z,y) 



The wave function is then deduced by solving the Weinberg equation (10) for the one boson exchange potential. In a 
weakly bound system, we can neglect to leading order the M 2 dependence of this potential by using M 2 w 4m 2 , and 
thus the expressions derived above for skewed parton distributions are valid. 

Next we must extend the definition of the potential by crossing quark lines. First we shall consider the case of one 
negative momentum fraction as required by the SPD. Within the Wick-Cutkosky model, we represent the interaction 
blob in Figure || by the exchange of one massless boson. We thus use the light-front time ordered graph to calculate 
the potential for quark pair production off the initial quark line. The interaction which describes hadronization off the 
quark line is then described by the graph in Figure |^.^ Notice we omit the diagram with pair creation from the vacuum. 

4 Notice that we do not include a spectator quark line of momentum q, say. This is because the time-ordering of the Z-graph is fixed by 
the q~ integration. Indeed we are merely interpreting the non-wave function vertex. Were we to include an intermediate state with energy 
q" in the denominator of Eq. B2, the crossed potential would have the wrong functional dependence. 
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Figure 8: Diagram of the crossed-pair potential (hadronization off the quark line) using one boson exchange in the 
Wick-Cutkosky model. The quark line at the top is altered to denote traveling backwards in time. 



The quark with momentum p, has y < and thus travels backwards in time (as the diagram indicates). Thus using 
the time-ordered rules, we arrive at Eq. ^2] for the crossed potential. Furthermore, we can use Eq. ^4| for the crossed 
potential since the contribution from Vf, (which would correspond to pair creation from the vacuum) vanishes. 

Next we need the potential for one momentum fraction greater than one. If the quark with z = p + / P + > 1, emits a 
4> boson which produces a quark pair, then one of that pair has the negative momentum fraction 1 — z and hence travels 
backwards through time. Using the time-ordered rules for this process yields an expression identical to 14 (where now 
V a vanishes due to vacuum pair production). Thus again we can use Eq. ^4] this time for for z > 1. It should not be 
surprising that the light front potential automatically contains such processes with "partons moving backwards in time" 
since the theta functions present in Eqs. (|52|, |53| ) do not restrict the momentum fractions to be between zero and one. 

Having the necessary potentials at hand, we now attempt to find an approximate solution to the Weinberg equation 
( |io| ) for the one boson, exchange potential (|54|) using the non-relativistic scheme. This method was first used by 
Karmanov f3~i|| to determine an approximate solution to the Wick-Cutkosky model, and so we summarize the procedure 



here. Following the approximation used in section 7A_, we can reduce the Weinberg equation to 

i 

- - M 2 



x{l 



^M)^-V(\,0; X M) J 2{2 %%_ z) ^) (56) 



V'O'k ) = T. 2 7Z 77^—, 7T 



(k^ 2 +(2x - l) 2 m 2 + x(l - x)m 2 a 2 ) 2 1 + |2a; - 1| ' 



in the limit m — > oo. The eigenvalue M 2 can then be solved for algebraically in the non-relativistic limit. We find 
M = 2m — \mo? . We then deduce the ground state wave function in this approximation by rewriting Eq. ( |56| ) 

(57) 

h x{L — x)m*a*)" x t 1^ -H 

where N is the normalization constant specified by Eq. ([131). 

This solution and the ensuing non-relativistic approximation scheme are valid to O(alna) from the weak binding 
limit and Q(—t/m ) from the non-relativistic approximation. When we look at the approximate non-valence SPD's in 
this scheme (Eq. |47| ), they are at 0{a 2 In a) relative to the valence region — the logarithmic dependence appears since 
the integral for C C (C, t) in Eq. ^ diverges when a — > 0. This is just as we might suspect in a weakly bound theory — 
pair creation and subsequent interaction are suppressed. Nonetheless, Eq. [l7| does give us a sense of the non- valence 
distributions, however, the sum rule Eq. can be verified to leading order using only the valence contribution to the 
SPD (Eq. pi| ). In Figure ||, we plot the percent contribution to the form factor from the valence SPD as a function of the 
skewness £. Lorentz invariance demands £ independence but our non-relativistic approximation retains Q dependence. 
The sum rule is satisfied by our non-relativistic approximation since the corrections are of 0(—t/m 2 ). 

Despite being suppressed by powers of a relative to the valence region, we can still look at the SPD's structure in 
the non- valence region (Eq. [l7]). The coefficient functions Cb(C,t) and C c (£,£) can both be calculated for small £ and 
—t/m 2 : we find Cb ~ +NQ 4 /% and C c ~ —N( 2 x 10 3 , where N — 0.00545 is the normalization constant appearing in 
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Figure 9: Verification of the sum rule in the non-relativistic limit. The percent error in the determination of the form 
factor (from the valence region SPD) plotted as a function of £ for small values of —t/m 2 . 




0.025 t/2 0.075 £ 



X 

Figure 10: Structure of the non-valence SPD in the non-relativistic limit —t/m 2 <C 1. Here we plot f(X,£) (which is 
proportional to J* JDw appearing in Eq. (^)) for the fixed value £ = 0.1. 



Eq. |5^ (determined with a = 0.08). Thus apart from variation with £ we can ascertain the non- valence SPD's behavior 
from / JDw In Figure [to] we plot f(X, Q = -± / ^j^I{<J, k x \t) D w (ct, kr 1 \t) as a function of X for fixed C- The 
choice of sign preserves the overall sign of T . The figure shows the SPD is symmetrical about X — £/2 arising from 
the symmetry of I and Dw about a — 1/2 (which in turn is due to the fact that our quarks are equally massive). Of 
course, as a function of a G (0, 1) the profile in the figure is universal. Given the nature of Cb and C c , since the scaled 
function f(X, £) is of order unity the non-valence SPD is ~ 10 6 times smaller than the valence SPD for small C and 
—t/m 2 (here we use £ = 0.1 for the comparison). The additional suppression (beyond powers of the weak coupling) is 
due to accessing the highly concentrated wave function near the end points. 

Using the non-relativistic limit to verify the sum rule is surely a good start, however, it does away with many of the 
complications inherent to the light-front approach. Were we to use the full solution to the Wick-Cutkosky model, not 
only would we have to worry about a consistent approximation scheme for the Green's function, but also the lack of 
covariance (here rotational invariance). The above expression Eq. (p2j) is exact and consequently covariant (although 
not manifestly so), but using any realistic procedure to solve for the wave function, we must truncate the interaction 
kernel. Doing so we loose covariance, which is approximately restored only when our approximations to the kernel are 
improved |3(|. This loss of rotational invariance would not only lead to improperly non-degenerate spectra but also 
would tamper with the sum rule. Without manifest Lorentz covariance, Eq. ( p2|) will only be approximately satisfied 
and a true challenge to verify exactly. This is a difficulty which plagues any Hamiltonian theory. 
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9 Conclusion and outlook 



Above we have seen how to express the pion's skewed parton distribution in terms of two body, light-front wave functions. 
By using a light-front reduced Bethe-Salpeter approach, we can express the SPD in terms of two body wave functions by 
crossing the interaction (and making the assumption that the effective qq potential is independent of energy). We did so 
first using only bare-coupling at the photon vertex. The functional form derived Eqs. (^4j27) is time reversal invariant 



(as we explained how to demonstrate through £ symmetry) and is continuous at the crossover. In this non-perturbative 



scheme, however, the bare-coupling photon vertex is not justified and thus we derived the photon vertex function Eq. 38 
in order to include the necessary qq interactions and hence the effect of higher Fock states. The interacting theory SPD 
( ff2| ) can then be written in terms of two body wave functions and the Green's function. The form derived maintains 
the necessary symmetry properties and is continuous. Potential schemes for approximating the Green's function were 
then investigated. We found the Born approximation to the photon vertex does not result in the bare-coupling SPD 
since the diagram corresponding to J- c Eq. (fid] ) is neglected. The non-relativistic approximation ( ffTj ) and the closure 
approximation (|5l]) to the non-valence region were then derived. 

Using the non-relativistic approximation, we were able to verify the sum rule for the Wick-Cutkosky model in the 
weak binding limit (see Figure^). Utilizing the crossed-pair potential and the non-relativistic approximation scheme, we 
explored SPD's in the non- valence region (see Figure |l^). This approach indeed provides an alternative to the explicit 
iV-body Fock space expansion, and hence a different interpretation of the contributions to the valence and non-valence 
regions. We believe this paradigm is the most useful intuitively, though currently suffers from calculational intractability. 
For instance: the mass squared dependence of the interaction has been neglected here (as is often the case elsewhere, 
e.g. constituent quark models) and remains a crucial issue to be dealt with, not only in the context of these distributions 
but for form factors and quark distribution functions as well. Additionally calculation of the full four-point Green's 
function presents another technical challenge. The closure approximation may however provide a way to reasonably 
approximate the Green's function, since continuity and the sum rule can be used as constraints to check how well the full 
Hamiltonian can be described by the two-body sector. Lastly, we find that use of the crossed-pair potential in tandem 
with the Green's function for the interacting theory allows one to obtain non- wave function vertices present in SPD's 
and time-like form factors. 
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Appendix A: £ symmetry 

At best, demonstrating the £ symmetry of our SPD's is subtle. To keep matters simple, we consider only the bare- 
coupling SPD's of section ^ and neglect quark self interactions, since their inclusion is straightforward and tangential to 
time-reversal invariance. While Eq. [2^ is clearly even in £, the sign of A + has been reversed and thus the location of the 
poles changed. Thus the whole derivation isn't necessarily valid for £ — > — £. Furthermore, directly switching £ — > — £ 
in equation ^9] yields zero. The time-reversal properties of our results are a crucial guide to consistency. Separately 
calculating the time-reversed diagrams of Figures | and ||, we easily arrive at the complex conjugates of Eqs. [29[ 
In order to show the relation between ^-symmetry and time-reversal, we recalculate the covariant triangle diagram with 
A-»-A. 

Looking at Figure ^ with A replaced with — A, we have the following expression for the covariant diagram 

d 4 k (2 k+ - A+)r(fc, P)T*(k - A, P - A) 

(2tt) 4 (fc 2 - m 2 + ie)((fc - A) 2 - m 2 + ie)((P - k) 2 - rn 2 + ie) ' (5 ' 

Having reversed the sign of A + in the diagram, we can again choose A + = — 2£P + < 0, where £ > and thus we have 
effectively taken £ here to be — £ of our previous calculation. Due to the kinematics, however, P + = (1 — £) P + . Thus 
defining A + = — £1 P + yields fi = This should be compared to the skewness £ = j+z- We still choose to work in the 
symmetrical frame, however, this constraint now reads P- 1 = A^/2 which is natural given the relative sign difference. 
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The poles of the k integrand are situated at 

= A" 



m 2 + (k- L -A^) 2 



2(k+-A+) 2(k+ -A + ) 

m 2 +k J ' 



2kT TP 



K = P~ - 



m 2 +(P^ -k 



2 



(P+-k+) ^ 2(P+ -k+) ' 



When k + < A + , all poles lie in the upper- half complex plane and when k + > P + , all poles lie in the lower-half plane. 
Thus, by Cauchy, these contributions to the integral vanish. Now we consider the non-zero contributions. 

Firstly we unravel the familiar case: < k + < A + which gives a contribution 2iri Res(fc~). At this point, however, 
the region appears far from familiar. The ratio k + / P + is restricted to (0, 1). To relate with our previous notation, we 
must define k + / P + = Y — j^f which restricts x € (£, 1). 

Using the definition established in equation we merely strip away the k + integral as well as the coupling to 
uncover 



^(x,c,t) 



dk 1 - 1 T(Y, k x -YA^/2) 

?.(9.ir)3 Y(Y 4- Q)(1 - Y) op4- p- m 2 +k^ 2 m 2 + (k^-A^/2) 



2P+{P- - A") 



m 2 + (k- L — A^) 2 _ mH(k x -A-L/2)^ 
Y+U l-Y 



(59) 



Simple algebra gives the relation = f+f = x - As before we can use A" = *+^ 2 to write 2 P+ A" = -(t+ A ±2 ) /CI. 
(Back in section [|, this relation was strikingly similar 2P + A - = —(t + A ±2 )/(.) We then exploit A • P = to find 



(1 + Cl/2) 2 A ±2 + Cl 2 M 
1 + Cl 



(60) 



which relates to the expression earlier by CI — > — (. Using the above, we can write 2P + A = Cl(M 2 + jA ±2 )/(l + CI). 
(Back in section g, we used the relation 2P+ A" = ((M 2 + |A ±2 )/(1 - £).) 
Identification of the initial meson's wave function in Eq. 59 is simple 

T(YM-YA-/2) _^ k x_ yA , /2) . (61) 



9 p+ p_ _ m 2 +k^ 2 _ m 2 + (k^ -A^/2) 2 
Zr r Y l-Y 



The final meson's wave function requires us to manipulate the energy denominator 

2p+(p- - A-> - + < k ; n A± > 2 - - 2 + - - + < k ;;;^- i)A± » 2 ). ««, 

Upon shifting the integral, we arrive at 

^ x -<'<> '/l- ym tf^ 1 + (i-*)AW,k-,. (63) 

The pre-factors along with the conversion factor 1 — £ 

1 + f2 1 - £ 2 

( 1- ^y(y + n)(i-y) = (x 2 -<e 2 )(i -x) ^ 

identically line up to our expression and thus we see why the naive replacement £ — ► — £ was justified. 

Lastly we consider the region A~< k + < 0. Given our identification of Y as k + /P + , we see x is restricted to 
( — £)£)• To evaluate the integral Eq. [5^ in this regime, we calculate — 2ni Res(fc~). From this expression, we toss away 
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the k + integral to uncover 



r(Y,k i -yA 1 /2) 



3 Y(Y + fi)(l - Y) 2 P + A- - "W" + m2+ ^~ A ") 2 



r*(Y,k^+(Y/2-l)A^) 



op+fp- _ a_s _ nf+^-A^Y _ m^ + (k^ -A-L/2) 



Comparing with our above equation (|59|) the final state meson vertex and energy denominator are unchanged and can 
thus be replaced by (1 + fl)if)*(X, k x +(X/2 - l)A x ). 
A little algebra needs to be employed to see 



r + y + n v ( y_ 1} 



1/ ™ 2 + (k ± +( < T-l)A- L ) ; ' 



where we used the relation Y/Vt = X/C, — 1 = a — 1. 

Now in Eq. [35], the initial meson's vertex can be related to a wave function only by using the interaction's crossing 
symmetry (since Y < 0) . Writing the remaining vertex in terms of the interaction, as well shifting the k^ integral by 
A ± we see 

-Q(l + fi) f d^dZd^ ± 
^ X '^> - Y(Y + n)(l-Y) J W^WW^z) A ' 



x V[Z, k ; Y, k +(1 - r/2)A- L J D w (a,k^ +aA ± \t)^*(X,k ± +XA ± /2). (67) 
For suspense, we have saved the pre-factors (and conversion factor) till the very end 

n n "(! + ") = ^(l+g) ffi8l 
1 ^>Y(Y + fi)(l-y) (z 2 -£ 2 )(l-:r)' 1 ; 

The £ symmetry of equation |2^ has been demonstrated and thus the covariant approach naturally incorporates the 
subtlety of switching the qq pair annihilation from the initial to final state when £ — ► — £. This algebra can easily be 
extended once interactions have been added at the photon vertex. 



Appendix B: SPD's, form factors and quark distributions 

Above we have made reference to verifying the sum rule Eq. (53) without specifying what structure the form factor 
has. Here we connect our expressions for SPD's to the celebrated Drell-Yan-West formula for the form factor |}7| and 
then analyze what neglecting the energy dependence of the qq interaction has done to SPD's, form factors and quark 
distributions. 

Since Lorentz covariance constrains the X integral of the SPD's to be independent of £, we may find the electromag- 
netic form factor by integrating the SPD's in the limit of zero skewness (£ — > 0). This corresponds to using the highly 
simplifying Drell-Yan-West frame in which the plus component of the photon's momentum vanishes (here A + = 0). This 
limit is now complicated by possible singularities arising from the non- valence region (X < £), like those encountered 
in j2TJ. Ignoring this complication for the moment, we see the contribution from the valence region (X > f) gives 

limj\r a (X,C,t)dX = J ^0^-^^xMW(xM+{l-X)A^) (69) 

the Drell-Yan-West formula. 
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It is interesting to first consider the non- valence contribution to the form factor coming from the bare-coupling SPD 
(Eq. |2^). As this region extends from X — to Q, to avoid singular momentum fractions, we change variables from X 
to a — X/C, and then take the limit in which £ vanishes. The contribution to the form factor reads 

^ . v . f da e?k^ dZ 'd k 1 - , ,. 

r^(Co-,k x ) 



limjf XT(X,(,t)dX ~ C J 



[2(2nmi-a)Z(l-Z) 

x A w (ff.^+dA 1 tWz./s" 1 ;^- l),k ± +A ± )^*(^,« ± ) (70) 



which vanishes as C ip(C er, . . . ) as £ — > 0. Thus for the Born approximation to the photon vertex, the Drell-Yan-West 
formula indeed gives the form factor. 

We can, of course, employ the same change of variables technique to consider the contributions to the electromagnetic 
form factor from the true SPD's of section |(| The contribution from Tk(X, £, t) Eq. ( |39| ) proceeds very similarly to the 
analysis for the bare-coupling SPD above. Additionally we must change variables from W £ (0, £) to u> = W/C, but in 
the end / X T h dX ~ ( V(C . . . ) as C 0. 



Lastly we consider the possible contribution from JF C (X, £, t) Eq. (41). The change of variables from X to £ runs 
as above. On the other hand, W £ (£, 1) giving us no reason for an addition change to ui. Taking £ to zero results in 

J X T c dX ~ £ y{^i . . . J where V is the crossed potential defined for a momentum fraction greater than one. 

Using the crossed potential for the Wick-Cutkosky model, this contribution to the form factor vanishes as £ 2 . 

We have shown that both contributions b and c vanish in the limit of zero skewness. Within the impulse approxi- 
mation, the Drell-Yan-West formula for the form factor is exact for our model. Thus we have seen, whether or not we 
include interactions at the photon vertex our model has the same form factor. The procedure here can be straightfor- 
wardly applied to the SPD's in the approximation schemes of section 0. For precisely the same reasons above, both 
schemes preserve the Drell-Yan-West formula in the forward limit. 

Having demonstrated there are no singularities associated with the limit C —> 0, we can also extract our model's 
quark distribution function q(x) = F(x, 0, 0). This is then simply the forward limit of Eq. ^J, namely 

<71) 

the square of the meson's wave function. In light of our above analysis, we know that the wave function vanishes 
quadratically at the end points and thus q(x) vanishes there as well. This is physical at x — 1 since one parton cannot 
posses all of the meson's plus momentum. At x — 0, however, the non-valence partons can rescue one parton from 
carrying all the plus momentum and so such configurations have non-zero probability, i.e. q(0) ^ 0. Indeed the fact 
that our model distributions vanish at x — is tied directly to neglecting the mass squared dependence of the effective 
interaction in the two-body subspace. The non- valence probability amplitude (V'nvlV'nv) is simply 

(Vw|</W> = -^2<V>2|F|V> 2 >, (72) 

where (V^IVIV^) is the projection of the interaction on the two-body subspace (simply denoted as V above). 

Accounting for the mass squared dependence of the effective interaction allows for non-valence contributions to the 
quark distribution functions. There are further consequences, however, given the unifying role of SPD's. The Drell-Yan- 
West formula, for example, will have corrections (of particular importance to the form factor's x- integrand at small x), 
and adding interactions at the photon vertex will then change the form factor. It has already been noticed |38) that 
model deuteron wave functions have better covariance properties than their form factors calculated from Drell-Yan-West 
expressions. Lastly the most general consequence is to affect the SPD's, particularly at the crossover, e.g. F(£, (,t) = 
seen here (as well as in [jLl[ p"2| ) stems directly from neglecting the mass squared dependence. The crossover point thus 
probes a convolution of non- valence states for which this approach may possibly be apt in describing, however, it would 
require a reworking of the photon vertex that does not make use of the potential model's Green's function. (As for 
the light-cone Fock space expansion, calculation of wave functions at small x is rather complicated but contains rich 
information about the relation between different Fock components 
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Figure 11: The triangle diagram for the time-like, pion form factor 



Appendix C: GDA's and the time-like form factor 

Time-like form factors have been often avoided in the light-front approach due to the unavoidable presence of non-wave 
function vertices. This complication is mirrored in the scheme of [^): the time-like pion form factor does not have a 
direct decomposition in terms of pion Fock component overlaps alone. Here we show how the above analysis can be 
applied to time-like form factors; that is, we obtain expressions in terms of two-body Bcthc-Salpeter wave functions and 
the four-point Green's functions. Furthermore, we make the connection with the generalized distribution amplitude for 
our model which encodes the non-perturbative physics of two pion production ^5) . We attempt below to preserve 
the notation set up in these references. 

The time-like pion form factor F(s) is defined by (see Figure |lTJ) 

MpMp') |J"( ) |0) = (p-p')"F(i), (73) 

where s = (p + p') 2 is the center of mass energy squared. Here the pion is merely our qq pair. Now define P M = p M + p' M 
and C = P + / P + - We can work out the kinematics of this reaction to find 

P S 



P = 



2P+ 

(i-C)« 



2P+ 

p X2 = s(l-CK-^ 2 , (74) 

where M is the pion mass and we work in the frame in which = 0. . 

The GDA for our model has a definition in terms of a non-diagonal matrix element of bi-local field operators 

*(*,C,«0 = J 27r(2z _ 1} (n(p)Ap')\q(x-)d + q(0)\0). (75) 



Such a definition of the GDA leads us immediately to the sum rule 

2z-l 



2C-1 



$(z,C,s)dz = F(s), (76) 



and hence a means to calculate $ from the integrand of the time-like form factor. 

Returning to Figure [H], we can write down the expression for the time-like form factor 

-i J -^T+(k,P)G (k)T(k,p)G (k-p)T(k-p,P-p)G (k-P) (77) 

and perform the light-front reduction by integrating out the minus dependence. Whenever troublesome momentum 
fractions appear, we use the equation of motion (110) to insert the crossed interaction. Lastly we merely extract the 
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Figure 12: Graphical representation of contributions to the GDA Eq. (|79[): $ = 4> a + 4>b- The shaded box again 
represents the crossed interaction and the internal double line represents the four-point Green's function. 



derivative coupling and z-integral from the photon vertex function to uncover the GDA. Using w = k + / P+, the relevant 
contribution from the photon vertex appears as 



(2z- 1)P+ 

with z £ (0,1). 

Carrying out this evaluation yields 



$(z,C,s) 



d q -1- dxud dxd k a 



2(2tt) 3 z(1 - z) 



-tp(x,K )G(z,q ;w,k |s) 



(78) 



[2(2n) 3 ] 3 w(w - C)(l - w)x(l - x)z{\ - z) 

x ( (f) a (w, k^; x, k X ) 0[w(£ — w)] + (f>b(w, k x ; x, k 1 ) (9[(1 — w)(w — £)] ] (79) 



where we have abbreviated 

Q («,, k-L ; x, k 1 - ) = C K (x,^ ; y , k x + (1 - y) p x ) ^ (| , k^ - 1 p J 

k x ; x, « x ) = (C - 1) V(x, k X ; k x — p X )^(y, k x +(1 - y) p x ) , 



(80) 



with y — j-^ ■ The two contributions (a and b) to the GDA can be interpreted graphically, see Figure [T^. 

We can furthermore decompose the GDA in the approximation schemes set forth in section fjj. In the non-relativistic 
scheme, we use the approximation to the Green's function in Eq. n& to derive 



$(z,C,s) 



NR 



dq A 



2(2tt) 3 z(1 - z) 



I(z,q X \s) D w (z,q x |s), 



(81) 



where 



^ r = 0) v(i,o^,-Eif W l,-Ei) 



K ,. ) ._^ I _„ ) v(l,^i.-^)*(i^.-^, 



(82) 



where i/)(r — 0) is defined in Eq. |4^. In this form, we see the non-relativistic approximation maintains continuity at 
( = 1/2. Apart from ( dependent coefficients, the structure of the non-relativistic GDA is identical to that of the 
non-relativistic SPD in the non- valence region (with a replaced with z), see Eq. (E^). 

Lastly appealing to the closure approximation for the Green's function, we arrive at the approximate GDA 



$(z,C,s) 



CL 



dq 1 - dxd K- 



-ipix,^) 



[2(2^) 3 ] 3 z(l-z)(z-CMl -x) 

I a (z,q- L ;£,Kr L ) 



0[*(C - z)\ + W *)(* C)]J ■ (83) 



22 



Notice given the form of the functions |8(J, the closure approximation to the GDA is continuous at z = C- This is a 
curious fact (in contrast to SPD's) since the factorization theorem for GDA's does not depend on continuity at z — ( 
p5| . Additionally at large s, we can apply the Born approximation to the Green's function. The resulting GDA's are 
identical to bare-coupling GDA's, which is a further contrast to SPD's. 
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